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ABSTRACT 

We report the Herschel SPIRE detection of dust emission arising from UV- luminous (L > L*) 
star-forming galaxies at 3.3 z ^ 4.3. Our sample of 1,913 Lyman Break Galaxy (LBG) candidates 
is selected over an area of 5.3 deg 2 in the Bootes Field of the NOAO Deep Wide-Field Survey. This 
is one of the largest samples of UV-luminous galaxies at this epoch and enables an investigation of 
the bright end of the galaxy luminosity function. We divide our sample into three luminosity bins 
and stack the Herschel SPIRE data to measure the average spectral energy distribution (SED) of 
LBGs at far-infrared (FIR) wavelengths. We find that these galaxies have average IR luminosities of 
(3 — 5) x 10 11 L Q and 60 — 70% of their star-formation obscured by dust. The FIR SEDs peak at 
A rcst £ lOO^im, suggesting dust temperatures (T d = 27 — 30 K) significantly colder than that of local 
galaxies of comparable IR luminosities. The observed IR-to-UV luminosity ratio (IRX = Lir/Ltjv) 
is low (« 3 - 4) compared with that observed for z « 2 LBGs (IRX Z ^2 ~ 7.1 ± 1.1). The correlation 
between the slope of the UV continuum and IRX for galaxies in the two lower luminosity bins suggests 
dust properties similar to those of local starburst galaxies. However, the galaxies in the highest 
luminosity bin appear to deviate from the local relation, suggesting that their dust properties may 
differ from those of their lower-luminosity and low-redshift counterparts. We speculate that the most 
UV luminous galaxies at this epoch are being observed in a short-lived and young evolutionary phase. 
Subject headings: galaxies: high-redshift — infrared: galaxies — ISM: dust, extinction 



1. INTRODUCTION 

Measurements of the spectral energy distribution 
(SED) of high-redshift galaxies can provide insight 
into their constituents (i.e., stars, gas, dust, AGN), 
star- formation rates (SFRs), and formation history. 
While high-redshift (z £ 2) star-forming galaxies have 
been identified in large numb ers using rest-frame UV 
color selection techniques fe.g.. iSteidel et al1 ll996. 2003; 
lAdelberger et al.ll2004t iBouwens et al.ll2010D . their faint- 
ness across the electromagnetic spectrum has limited our 
understanding of their multiwavelength properties. In 
particular, UV emission is easily extinguished by dust 
and re-radiated at far-infrared (FIR) wavelengths. With- 
out direct measurements of the FIR emission, the amount 
of dust extinction and the intrinsic SFRs remain uncer- 
tain. Thus far, our understanding of the SFRs largely 
rests on the assumption that the UV-selected high-z 
galaxies are similar in their dust properties to local star- 
burst galaxies, and specifically that the dust extinction 
in the UV (as measured by the slope, (3, of the UV con- 
tinuum) can be used to predict th e infrared luminos- 
ity (Lir) due to dust emission (e.g., iMeurer et all 119991: 
Calz etti et alll200Ct iReddv et al 2008; Finkelstein et all 



2009). 
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The Spitzer and Hersche^\ Space Telescopes have en- 
abled direct measurements of the mid- and FIR emis- 
sion from distant galaxies and thus a direct test of 
this assumption. While these telescopes still lack the 
sensitivity to individually detect any but the most IR- 
luminous high-z galaxies, stacking analyses of z ~ 2 — 3 
UV-selec ted galaxy samples have r esulted in use f ul con- 
straints fe ifiopoulOTreLall 120101 : iMagdis et all I2010bt 
IReddv et all 120121) . Several studies of galaxies at z ~ 
2-3 ( Reddv et al.ll2006tlDaddi et a l.H2007t IMa gdis et all 
I2010at IReddv et alll2010l 12012ft have demonstrated that 
the ratio of the average UV continuum slope to the av- 
erage far-infrared luminosity is similar to that in local 
starburst galaxies. There are exceptions: most notably, 
there is evidence that the local /3-Ltr relation may not 
apply to very dusty, UV-faint galaxies, or to extremely 
youn g galaxies (e.g., iGoldader et all 120021 : IReddv et all 
120101) . Also, the studies thus far are based primarily on 
galaxy samples with luminosity ranges that do not in- 
clude large numbers of the most luminous galaxies, and 
the dependence of the /3-Lir relation on luminosity is not 
well understood at high-z. 

In this paper, we extend this test to higher-redshift 
and higher-luminosity UV-selected galaxies. We inves- 
tigate the far-IR properties of a very large sample of 
3 z ?i 4.5 galaxies. We stack the Hers chel Spectral 
and P hotometric Imaging Receiver (SPIRE: iGriffin et all 
120101 ) observations of this sample and present the first 
(stacked) detections of the far-IR emission of UV-selected 
galaxies at this epoch. We use a ACDM cosmology with 
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(n,fi A , 0-8,^100) = (0.28,0.72,0.9,0.72). 

2. DATA AND GALAXY SAMPLE 

The datasets used in our analyses are: optical 
data (B W RI) in t he Bootes field of the NOAO Deep 
Wide-Field Survey (jj annuzi fc Devl H999'): near-infrared 
JHK S data taken with the NEWFIRM camera (Gon- 
zalez et al., in prep); mid- and far-infrared (3.6- 
160/im) data from the Spitzer Infrared Array Cam- 
era (IRAC) and the Multiband Imaging Photometer 
(MI PS) obtained by t he Spitzer Deep Wide-Field Sur- 
vey (jAshbv et al.|[2009h and the MIPS AGN and Galaxy 
Evolution Survey (Jannuzi et al., in prep); and far- 
IR (250-500/im) Herschel SPIRE observations obtained 
by th e Herschel Multi-ti ered Extragalactic Survey (Her- 
MES: iQliver et al.ll2012h . We reduced and mosaiced the 
SPIRE data using t he Hersche l Interactive Processing 
Environment (HIPE; EH |2010j) , removing striping, as- 
trometry offsets, and glitches missed by the standard 
pipeline (Alberts et al., in prep). The SPIRE maps in- 
clude a deep 2 deg 2 region surrounded by a shallower 
outer region. The lcr depths at 250, 350 and 500/im are 
3, 2, 3 mjy in the inner region and 5, 4, 5 mJy in the 
outer region respectively. 

We use broad-band (By/RI) color selection to iden- 
tify sources that e xhibit a continu um break between the 
By/ and R bands dLee et al.ll2011[) This technique, pio- 
neered bv lSteidel fc Hamiltori l 19931) . successfully identi- 
fies z ~ 4 UV-emitting star-forming galaxies (also known 
as Lyman Break Galaxies, or LBGs) which exhibit a 
strong break in their rest-frame UV spectra at 912A (the 
Lyman limit) and additional flux attenuation due to the 
intervening Lyman alpha forest at 912 — 1216A. Spectro- 
scopic redshifts of ~ 5% of the sample (|Lee et al.l 120111 : 
iKochanek et al.ll2012L Lee et al. 2012, in prep) have con- 
firmed a broad redshift distribution (3^2^, 4.5, z m 
3.7; <7 2 s=s 0.4), and a contamination rate of 8 — 23% 
(the limits representing the optimistic/pessimistic sce- 
nario based on the nature of the ambiguous spectra). 
Our final sample consists of 1,913 Z/tjv > L* LBG can- 
didates selected over a contiguou s 5.3 deg 2 region of the 
Bootes field fsee lLee et al~ll2011l ). 

We divide the sample into three bins ac- 
cording to the /-band magnitudes as Jab = 
[21.3, 23.7], [23.7, 24.3], [24.3, 24.7]. For reference, 
the characteristi c luminosity L* at z = 3.7 corresponds 
to Jab = 24.7 (|Bouwens et al.l 120071 ). The magnitude 
bins correspond to overlapping luminosity ranges due to 
the finite width of the redshift distribution. By using 
the observed magnitude and redshift distributions of the 
sample, we find that the magnitude bins correspond to 
roughly Gaussian bins of Ltjv / L* , with widths at half- 
maxima of [>2.3], [1.4,2.5], and [0.9,1.5], respectively. 
We refer to these bins as the "high-", "intermediate-", 
and "low-luminosity" sample. The majority (97%) of 
galaxies in the high-luminosity bin are / > 22. The 
median UV luminosities of our samples are presented in 
Table [U 

3. DATA STACKING 

We measured the average optical-to-mid-IR SED of the 
galaxy samples by using photometry of stacked images in 
the individual optical, near-IR and IRAC bands. We first 



convolved all images in a given band to a common point- 
spread function (PSF), then extracted a cutout centered 
on each source, subtracted a local sky background, and 
produced a composite image by taking the median flux 
ove r all sources at e ach pixel position (for further details, 
see lLee et al"]l2011[ ). 

The Spitzer MIPS and Herschel data both suffer from 
large beams and spatial variations in the sensitivity. For 
the MIPS data, we used a variance-weighted median to 
account for small-scale variations in the S /N of the im- 
ages as a function of position. The stacked data only 
yielded upper limits. The 3cr limits, measured in aper- 
tures of diameter 6", 18", 38" (matched to the full-width- 
at-half-maximum of the PSFs in the 24^m, 70fj,m and 
160/im images, respectively) are presented in Table [2j 

The HerMES SPIRE data only cover a portion of our 
field, and sample 1673 LBGs (i.e., 85%). Thirty five 
(2%) LBGs have a 250/zm detection (> 5a) within 8"; 
33 of these have more than one optical counterpart in 
the beam. Of these, only two sources exhibit increasing 
flux densities towards longer wavelengths, as would be 
expected for z > 3.5 galaxies (sampling the Wien side 
of the dust emission). However, both have another op- 
tical counterpart other than the LBG candidate within 
the beam, making it uncertain whether the far-IR emis- 
sion is uncontaminated. We therefore assumed that the 
detected sources are contaminated by unresolved inter- 
lopers, and excluded all 35 sources from the stacking. 
We stacked the SPIRE data by calculating the variance- 
weighted average flux. The measurement significance 
was d etermined using a bootstrap resampling technique 
(e.g., iPress et all Il992). We measured the flux using 
10,000 random samples of the source list (replacing cho- 
sen sources each time) and measured the RMS from the 
resulting flux distribution. Our results are presented in 
Tabled For completeness, we repeated our analyses with 
and without the 250/j,to detected sources; the differences 
are within the lcr uncertainties. 

High-redshift submill imeter galaxies are known to 
be st rongly clustered (jBlain et alJ 120041 : iHickox et al.1 
2012) and are often fo und in close proximity to LBGs 
([Chapman et al.l 120091 ). The large SPIRE beam could 
be contaminated by such sources. In order to estimate 
the degree of contamination, we used the substantially 
deeper 24^tm MIPS data and measured photometry in 
large apertures matched to the SPIRE beamsize. We 
carried out this measurement at both on-source (LBGs) 
and random sky positions. Wc found that (1) there is no 
significant difference between the on-source and sky dis- 
tributions, and (2) the distributions are similar at both 
the 250/im and 500^m beamsizes. This suggests that the 
contribution of clustered sources to the flux measured for 
the LBGs is small. Therefore, our stacked fluxes likely 
represent the FIR emission from a typical z ~ 3.7 LBG. 

4. RESULTS AND DISCUSSION 

4.1. The Panchromatic View of U\ 7 -Luminous 
Star- Forming Galaxies at z ~ 3.7 

Figure Q] shows the stacked SEDs of the LBGs in the 
three magnitude bins. At A ros t < Wfj,m, the SEDs show 
similar shapes, with two noticeable differences. First, the 
UV continuum slope /3 (f\ oc A^ measured at A rcst = 1300- 
3000A) increases with luminosity from /3 = —1.75 ±0.15 
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TABLE 1 

Summary of UV-to-FIR properties of Star-Forming Galaxies atz~ 3.7 



Sample N gal a log[L uv ] b log[Li R ] b IRX SFR UV C SFR IR C 

7ab = [21.3.23.7] 325/293 11.21 ±0.06 (11.65 - 11.74) 2.7- 3.3 -1.24 ±0.08 49 ± 7 77-95 

i" A B=[23.7,24.3] 784/686 10.96 ± 0.06 (11.48 - 11.57) 3.4-4.1 -1.65 ±0.11 28 ± 4 52 -64 

i" A B=[24.3,25.0] 794/694 10.76 ± 0.06 < (11.50 - 11.57; 3<r) < (5.5 - 6.4) -1.75 ±0.15 18 ± 3 < (54 - 63) 



a The number of LBGs and the number of LBGs included in the Herschel stacking 
b In units of Lq 

c In units of AiQyr -1 . The Kcnnicutt (1998) calibration is used assuming the Salpeter IMF (0.1 — IOOMq), solar metallicity 



TABLE 2 

Flux Densities in the Stacked MIPS and Herschel Images 





Sample 1 


Sample 2 


Sample 3 




I AB =[21.3,23.7] 


7ab=[23.7,24.3] 


/ AB =[24.3,24.7] 


24fim 


< 0.016 


< 0.011 


< 0.011 


70/xm 


< 1.38 


< 0.90 


< 0.90 


160/im 


< 13.2 


< 9.54 


< 9.92 


250/im 


< 1.05 


< 0.84 


< 0.72 




< 1.11 


< 0.90 


< 0.78 


350/mti 


0.65 ±0.39 


0.44 ±0.29 


< 0.78 




1.01 ±0.43 


0.85 ±0.33 


< 0.87 


500/mti 


1.45 ±0.47 


1.06 ±0.31 


< 0.90 




1.69 ±0.48 


1.37 ±0.33 


< 0.96 



Note. — The flux densities are in units of mjy and the upper 
limits are 3cr in flux. For the Herschel data, the first and second row 
at each observing wavelength represent the stacked fluxes excluding 
and including sources with a > 5<r detection at 250fim, respectively 



(low-luminosity bin) to —1.24 ± 0.08 (high-luminosity 
bin). Shifting the median redshift changes the level of 
intergalactic attenuation and but does not impact the 
observed colors significantly (e.g., Az = 0.1 would re- 
sult in A/3 ;$ 0.05). Second, the break observed between 
J and K$ decreases in strength with increasing lumi- 
nosity. If this is the Balmer break, the trend suggests 
that galaxies in the high-luminosity bin are, on average, 
younger than their lowe r luminosity counterparts (e.g., 
IBruzual fc CharlotlH99l . 

In the far-IR, only the intermediate- and high- 
luminosity bins show detections in the 350/ito and 500/im 
SPIRE bands. To measure the total IR luminosities, we 
model ed the SEDs using the templates of lCharv"fe r Elbaz 
(|2001l hereafter CE01). The CE01 templates are de- 
signed to match the SEDs of present-day IR-luminous 
galaxies and to reproduce the observed local correlation 
between dust temperature (Td us t) and IR luminosity. If 
we only use CE01 templates that correspond to the ob- 
served IR luminosity ranges, the resulting fits are poor 
(Xr ^ 2.3; 90% confidence). This is primarily because 
the observed SEDs of our samples peak at longer wave- 
lengths than those of local galaxies of comparable IR 
luminosities. The longer peak wavelength observed for 
our galaxies may suggest that galaxies at z ~ 3.7 have 
colder dust than local galaxies at similar IR luminosi- 
ties. Similar results have been found for various types 
of high-redshift (z > 1 .5) galaxies ([Pope et al.l 120061: 
Dannerbauer et al.l[2fjlfl:lMuzzin et aLlboiOHElbaz et ail 
20111 : IKirkpatrick et al.ll2012t ). 

To account for possible redshift evolution of IR SEDs, 
we fit all the CE01 templates leaving the normalization as 
a free parameter, effectively eliminating the local Tdust- 
Llr correlation from the fitting procedure. The 3cr upper 
limits were included in the % 2 calculation only when the 



model exceeds the limit. For the low-luminosity bin with 
no Herschel detection, we determined the upper limit on 
the IR luminosity using the best-fit templates of the two 
other bins. The total IR luminosities (Lir: integrated 
between 8 — 1000/im) are computed by integrating the 
scaled best-fit CE01 template. Typically, \r ~ 0-1 f° r 
the best-fit models. 

We find L m values of (4.5 - 5.5) x 1O U Z , (3.0 - 3.7) x 
lO n L , and < (3.2 - 3.7) x lO n L , in the order of de- 
creasing UV luminosities (Table [lj . Despite the uncer- 
tainties, a clear trend emerges: more UV-luminous galax- 
ies are also more IR-luminous. The main source of uncer- 
tainty is the peak location of dust emission (and therefore 
Td). The lower limit corresponds to Td « 30K (mod- 
eled as a grey-body function /„ oc B v {Td,v)vP where 
B v is the Planck function, (9=1.5 is the dust emissiv- 
ity). The upper limit is set by the coldest CE01 template 
{T d « 21K). 

As all SPIRE bands lie at shorter wavelengths than 
the FIR peak at z ~ 3.7, we are unable to constrain 
the SED shape from our data alone. However, it is in- 
structive to consider the existing measurements at longer 
wavelength s on other UV - select ed LBG samples at z > 3. 
At z ~ 3 iWebb et all (|2003h reported a 1.5a detec- 
tion at 850 nm (using SCUBA) for L > L* LBGs, while 
iMagdis et alj (|2010aD found a 3.7c detection at 1.1mm 
(using AzTEC) of IRAC-detected LBGs. At z ~ 5, 
iDavies e~a l. (2012) determined a 2a upper limit of 1.6 
mjy at 1.2mm (using MAMBO). In Figure [U we show 
these points after correcting the wavelengths and flux 
densities to z = 3.7. Even though the different selec- 
tion methods and redshift ranges of these studies make 
direct comparisons difficult, the measurements still pro- 
vide a useful constraint on the SED, preferring a peak 
wavelength near « 500/Ltm. The inferred dust temper- 
ature is Td = (27 — 30) K. When compared to recent 
Herschel-b&sed templates, the combined FIR constraints 
are in g ood agreement with the z ~ 2 star- forming tem- 
plate of iKirkpatrick et al.l (j2012f) and the starburst tem- 
plate of lElbaz et al.l ()2011[ )~ The main-sequence Elbaz 
template is ruled out (> 2.5cr) as it predicts 3 — 4x 
higher flux densities in t he sub-mm/mm wavelengths (see 
IKirkpatrick et af]l2012L for further discussion). Further 
constraints can be easily obtained by ALMA. Our best- 
fit model for the 1.5L* - 2.5L* LBGs predicts the flux 
densities of 0.83 (0.51) mjy in the 450 (870)Aim, which, 
according to the ALMA Sensitivity Calculator, can be 
detected (5a) in 2.1 hr (5.6 min). 

4.2. Dust-Obscured Star Formation 

Having measured the IR luminosity of the LBG sam- 
ples, we now compare the relative energy budget of star- 
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Fig. 1. — Lef t: The median SEDs ( 0.3 — 500fim) of the three samples are shown with the best-fit stellar population models and local 
IR templates of Chary & Elbaz (2001). The Herschel points are slightly offset (in wavelength) for clarity; the intermediate-luminosity 
bin is positioned at the correct wavelengths. The peak of the dust emission is not well constrained by the current data and we therefore 
show the acceptable models with both the lowest /highest dust temperatures (Td = 21/30K). The latter model is favored by the sub- 
mm/mm measurements of similarly-selected LBGs (Webb ct al. 2003; Magdis ct al. 2010a; Davics ct al. 2012). Right: The same far-infrared 
measurements are compared with three Herschel-b&sed templates of Elbaz ct al. (2011) and Kirkpatrick ct al. (2013). All templates are 
normalized to match the intermediate-luminosity points (filled triangles). 

formation in the UV and infrared. For each bin, we esti- 
mate the luminosity ratio IRX= Ltr/Ltjv assuming two 
dust temperatures, T d = 30K and 21K (see Table [1]). 
In all luminosity bins, the IRX remains relatively low at 
IRX w 3 — 4, implying 60 — 70% of the star formation 
in these LBGs is obscured by dust. The IRX is weake r 
than that observed at z ~ 2 f»80%: iReddv et al.|[2012h , 
suggesting the buildup of dust with cosmic time as low- 
mass stars continue to evolve into the AGB phase (e.g., 
iBouwens et~alll2009| ). 

Next, we investigate how dust extinction inferred 
from the UV comp ares with dust emission in the IR. 
iMeurer et al.l (|1999D derived the scaling law between 
dust extinction and the UV slope /3 based on the ob- 
served correlation for local starbursts. Such a corre- 
lation is ex plained if local starb urst galaxies are gov- 
erned by the lCalzetti et all (|2000f ) extinction law, which 
gives a natter attenuation curve than that found for the 
Large and Small M agellanic Clouds (jFitzpatrickl 119861 : 
iBouchet et ahlll985l). 

At z ~ 2. iReddv et all (|2010L 120121 ) found that most 
UV-selected star-forming galaxies follow the local IRX-/3 
relation, demonstrating that the dust properties and rel- 
ative geometry between dust and stars within the z <~ 2 
star-forming gal axies are similar to those in loca l star- 
bursts (see also iNandra et al.l [200l iDaddi et all 120071 : 
iMagdis et al.ll2010al ). However, there are also exceptions: 
several galaxies lie off this relation (Reddy et al 2006, 
2012; Siana et al. 2008, 2009) and are interpreted as 
being very young (<100 Myr). 

Figure [2] shows the locations of the three z ~ 3.7 LBG 
subsamples on the I RX-/3 plane and the local relation of 
IMeurer et all (IT999f1. For comparison, we also show the 
locatio ns of the z ~ 2 L\jy « L* galaxies of IReddv et al.l 
(2012) and t wo lensed galaxies (the Cosmic Eye and 
MS1512-cB58: ISiana et all [20081 [20091) . The galaxies in 

7 The correlation is scaled by 0.24 dex to match our definition 
of Ltr as IR luminosity integrated at 8 — lOOOfim from that of 



Mcurcr ct al. (1999, at 40 - 120/im) 



Fig. 2. — The IRX values as a function of the UV slope /3: The 
local Meurer relation, is shown with the 0.4 dex scatter. We show 
the locations of three subsamples at z ~ 3.7 assuming two different 
dust temperatures (T^ = 30K, 20K in filled and open symbols, 
respectively) . Both estimates are slightly offset in f} for clarity. The 
locations of the z ~ 2 galaxies, young (< 100 Myr) z ~ 2 galaxies 
(IReddv et al.120121) . and two lensed LBGs (ISiana et al.l2008l.l2009h 
are also shown. The galaxies in the two less-luminous bins follow 
the local relation while the galaxies in the high-luminosity bin lies 
outside the nominal range the Meurer relation close to young z ~ 2 
galaxies. 

the intermediate-luminosity bin closely trace the local re- 
lation. The upper limits for our low-luminosity sample 
are also consistent with the local relation. Hence, our 
results suggest that the dust properties in the L* — 2.5L* 
LBGs at z ~ 3.7 are similar to those of local starbursts. 

The high- luminosity bin lies significantly (0.5 dex) be- 
low the Meurer relat ion, closer to the location of the 
young z ~ 2 galaxies (jReddv et al.ll2012h . The observed 
offset suggests that the dust properties of the galaxies 
in the high-luminosity bin may be different from those 
of their lower-luminosity counterparts. Also, since the j3 
values for these galaxies cannot predict the LrR using the 
local IRX-/3 relation, it implies that SFRs derived from 
the UV using the standard relationship overpredict the 
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true SFR in these systems. The impact of this result on 
the total cosmic SFR density is minimal. Nevertheless, 
the offset in the IRX-/3 plane has interesting implications 
for the evolutionary state of these high-luminosity galax- 
ies. 

The high-luminosity LBGs appear to have a weaker 
Balmer break than those in less-luminous bins, im- 
plying that they may h ave younger population ages. 
iReddv et al.l (|2010t |2012|) made similar observations that 
z ~ 2 galaxies with ages < 100 Myr have lower IRX val- 
ues. This offset suggests that the most UV-luminous 
galaxies may be undergoing a short evolutionary phase 
during which the dust covering fraction and/or the dust 
grain size distribution is different than that in their less 
luminous counterparts (Siana et al. 2009, Reddy et al. 
2010). 

The statistics from the existing spectroscopy supports 
the possibility that a large dust covering fraction may be 
suppressing the emergence of Lya photons in the most 
UV-luminous galaxies. Of the 13 galaxies at Iab = 
22.5 — 24.3, none has strong Lya emission (rest-frame 
equivalent width EWq <20A). In contrast, 12/30 galax- 
ies at Jab = 24.3 - 25.0 have EW > 20A. While the 
spectroscopic confirmation of the fainter LBGs may be 
biased toward the stronger Lya emitters, the lack of Lya 
emission among the most UV-luminous galaxie s is signif- 
icant (see IStark et al]l20l3 IKornei et al.l[2010l ). 

The suggested luminosity-age relation we observe in 
the z ~ 3.7 LBGs is not observed for z ~ 2 galaxies 
as the median luminosity of their youngest (IRX-outlicr) 
population is similar to that of the overall sample. How- 



ever, the luminosity range of our high-luminosity sample 
is not well represented in the Reddy et al. sample. It is 
also possible that our results may have been influenced 
by sample contamination. Interlopers such as dusty low- 
redshift (z < 1) galaxies or low-luminosity quasars could 
simulate a smaller Balmer break and depressed FIR flux. 
Based on our DEIMOS spectroscopy, the contamination 
rate in the high-luminosity bin is «15%, slightly higher 
than that of the full sample (10-12%). Even if the FIR 
flux is zero for the interlopers, the resulting UV flux cor- 
rection to the IRX would be insufficient to explain the 
observed offset in Figure 2. If the low-redshift interlopers 
contribute disproportionately to the measured average 
IR flux, the true IRX offset of our high-luminosity LBGs 
from the Meurer relation would be even larger than ob- 
served. Spectroscopic confirmation of a larger number 
of candidates can place more stringent limits on the in- 
terloper populations, and also allow the stacking of con- 
firmed LBGs. 
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